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ABSTRACT 



Aims. We reexamine the theoretical instability domain of pulsating PG1159 stars (GW Vir variables). 

Methods. We performed an extensive g-mode stability analysis on PG1159 evolutionary models with stellar masses ranging 
i from 0.530 to O.741M0, for which the complete evolutionary stages of their progenitors from the ZAMS, through the thermally 

■ pulsing AGB and born-again phases to the domain of the PG1159 stars have been considered. 

Results. We found that pulsations in PG1159 stars are excited by the K-mechanism due to partial ionization of carbon and 
oxygen, and that no composition gradients are needed between the surface layers and the driving region, much in agreement 

■ with previous studies. We show, for the first time, the existence of a red edge of the instability strip at high luminosities. We 
"^h ' found that all of the GW Vir stars lay within our theoretical instability strip. Our results suggest a qualitative good agreement 

C*) , between the observed and the predicted ranges of unstable periods of individual stars. Finally, we found that generally the 
I 1 seismic masses (derived from the period spacing) of GW Vir stars are somewhat different from the masses suggested by 
O | evolutionary tracks coupled with spectroscopy. Improvements in the evolution during the thermally pulsing AGB phase and/or 
+3 ■ during the core helium burning stage and early AGB could help to alleviate the persisting discrepancies. 

ctf ■ 

^ Key words, stars: evolution — stars: interiors — stars: oscillations — stars: variables: other (GW Virginis) — white dwarfs 

h : 

. 5r 1 !■ Background GW Vir stars (those without nebulae) 1 . GW Vir stars 

are particularly important to infer fundamental properties 

Pulsating PG1 159 stars — after the prototype of the spec- about pre - white dwarfs in general, such as the stellar mass 

tral class and the variable type, PG 1159-035 or GW Vir and the sur f ace compositional stratification 2 (Kawaler & 

- are very hot hydrogen-deficient post- Asymptotic Giant Bradley 1994; Corsico & Althaus 2006). 
Branch (AGB) stars with surface layers rich in helium n/^ncn i. uvji-Ui-i. 1 ±- 

30 85/0 b ( 15 60^) d 2 PC1159 stars are believed to be the evolutionary 

) oj, car on|~ oj an oxygen [~ °) connection between post-AGB stars and most of the 

(Werner & Herwig 2006) that exhibit multiperiodic, low , , , c ■ , 1 r ™ , , , 

\ , . ° ' * ' hydrogen-deficient white dwarfs. Ihese stars are thought 

degree [l < 2), high radial order [k > 18) o-mode lumi- , , ii.ru • • j A ■ j 

° v — , ' ' . i > to be the result of a born again episode triggered either 

nosity variations with periods m the range from about 300 by & yery kte helium thermal pulge (VLTp) occurring 

to 3000 seconds. Some GW Vir stars are still embedded - n & hot white dw&rf ghortly &fter hydrogen burning has 

in a planetary nebula; they are commonly called PNNVs almogt ce&sed (gee Fujimoto mi . Sch o nbe rner 1979 and 
(Planetary Nebula Nucleus Variable). PNNV stars are 
characterized by much higher luminosity than the "naked" 



In this paper, naked GW Vir stars and PNNV stars will be 

indistinctly referred to as GW Vir stars or pulsating PG1159. 

Send offprint requests to: A. H. Corsico 2 In addition, pulsating PG1159 stars have recently been 

* Member of the Carrera del Investigador Cientffico y shown by Corsico & Althaus (2005) to be valuable tools to 
Tecnologico, CONICET, Argentina. constrain the occurrence of extra mixing episodes in their pro- 

** Fellow of CONICET, Argentina. genitor stars. 
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more recently Althaus et al. 2005) or a late helium thermal 
pulse (LTP) that takes place during the post-AGB evolu- 
tion when hydrogen burning is still active (see Blocker 
2001 for references). During a VLTP episode, most of the 
hydrogen-rich envelope of the star is burnt in the helium- 
flash convection zone, whilst in a LTP hydrogen-deficient 
composition is the result of a dilution episode. In both 
cases, the star returns rapidly back to the AGB and fi- 
nally into the domain of high effective temperature as a 
hydrogen-deficient, quiescent helium-burning object. 

A longstanding problem associated with pulsating 
PG1159 stars is related to the excitation mechanism. The 
early work by Starrfield et al. (1983) was successful in 
finding the correct destabilizing agent, namely the re- 
mechanism associated with the partial ionization of the 
K-shell electrons of carbon and/or oxygen in the enve- 
lope of models. However, their models required a driving 
region very poor in helium in order to be capable to ex- 
cite pulsations; even very low amounts of helium could 
weaken or completely remove the destabilizing effect of 
carbon and oxygen (i.e. "helium poisoning"). The latter 
requirement led to the conjecture that a composition gra- 
dient would exist to make compatible the helium-devoid 
driving regions and the helium-rich photospheric compo- 
sition. Even modern detailed calculations still point out 
the necessity of a compositional gradient in the envelopes 
of models (Bradley & Dziembowski 1996; Cox 2003). The 
presence of a chemical composition gradient is difficult 
to explain in view of the fact that PG1159 stars are still 
experiencing mass loss (M ~ 10~ 81 A/q yr -1 for PG 1159- 
035; Koesterke et al. 1998), a fact that prevents the ac- 
tion of gravitational settling of carbon and oxygen, and 
instead, tends to homogenize the envelope of hot white 
dwarfs (Unglaub & Bues 2000). 

Clearly at odds with the hypothesis of a composition 
gradient in the PG1159 envelopes, calculations by Saio 
(1996), Gautschy (1997), and Quirion et al. (2004) - 
based on modern opacity OPAL data — demonstrated 
that g-mode pulsations in the correct ranges of effec- 
tive temperatures and periods could be easily excited in 
PG1159 models having an uniform envelope composition. 
The most recent study about PG1159-type pulsations is 
that of Gautschy et al. (2005) (hereinafter GAS05) based 
on a full PG1159 evolutionary sequence started from the 
zero-age main sequence (ZAMS) and evolved through the 
thermally pulsing and VLTP phases (see Althaus et al. 
2005). These authors found no need for invoking com- 
position gradients in the PG1159 envelopes to promote 
instability. 

As important as they are, the vast majority of the 
studies of pulsation driving in PG1159 stars rely on 
simplified stellar models. Indeed, the earliest works em- 
ployed static envelope models and old opacity data. Even 
more modern works, although based on updated opacity 
data (OPAL), still use a series of static envelope models 
that no represent a real evolutionary sequence, or evo- 
lutionary computations based on simplified descriptions 
of the evolution of their progenitors. The only excep- 



tion is the work of GAS05, which employs equilibrium 
PG1159 models that evolved through the AGB and born- 
again stages, beginning from a 2.7M Q zero age main se- 
quence model star. GAS05 analyzed four model sequences, 
with 0.530,0.55,0.589 and O.64M , being the 0.589M Q 
sequence derived directly from the evolutionary computa- 
tions of Althaus et al. (2005). The remainder sequences 
were created from the 0.589M Q one by appropriately 
changing the stellar mass shortly after the end of the born- 
again episode. 

On the basis of full evolutionary PG1159 models cov- 
ering the whole range of observed GW Vir masses, this 
paper is intended to confirm and extend the results al- 
ready put forward by the stability analysis by Saio (1996), 
Gautschy (1997), Quirion et al. (2004), and GAS05. We 
analyze the pulsational stability of seven different evolu- 
tionary sequences of PG1159 models with stellar masses 
between 0.530 and 0.741M©. Here, all of the PG1159 evo- 
lutionary sequences have been derived by considering the 
complete evolution of their progenitors, an aspect that 
constitutes an improvement over previous studies. One 
of such sequences (0.589M Q ) is that already presented 
by Althaus et al. (2005) and analyzed by GAS05, and 
the remaining ones are those computed recently by Miller 
Bcrtolami & Althaus (2006), with the exception of the 
O.741M0 sequence, which is presented for the first time 
in this work. We believe that the pulsational results pre- 
sented here based on extensive full evolutionary models 
shed new lights on the GW Vir stars and place previous 
studies on a solid basis, regarding stellar modeling. The 
paper is organized as follow: in the next Section we briefly 
describe the input physics of our evolutionary code and 
the PG1159 evolutionary sequences analyzed. A brief de- 
scription of our nonadiabatic pulsation code is presented 
as well. In SJ2I we elaborate on a detailed description of 
our stability analysis and in we compare our predic- 
tions with the observed properties of known GW Vir stars. 
Finally, in Sect. [5] we summarize our main results and 
make some concluding remarks. 

2. Evolutionary sequences 

The nonadiabatic pulsational analysis presented in this 
work relies on stellar models that take into account the 
complete evolution of the PG1159 progenitor stars. The 
evolution of such models has been computed with the 
LPCODE evolutionary code, which is described in Althaus 
et al. (2005). LPCODE uses OPAL radiative opacities 
(including carbon- and oxygen-rich mixtures) from the 
compilation of Iglesias & Rogers (1996), complemented at 
the low-temperature regime with the molecular opacities 
of Alexander & Ferguson (1994) (with solar metallicity) . 
Chemical changes are performed via a time-dependent 
scheme that simultaneously treats nuclear evolution and 
mixing processes due to convection, salt finger and over- 
shooting. Convective overshooting is treated as an expo- 
nentially decaying diffusive process above and below any 
convective region. 
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Table 1. Initial and final stellar mass (in solar units), and 
the final surface chemical abundances by mass (PG1159 
regime) for the evolutionary sequences considered in this 
work. 
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0.530 


0.33 


0.39 


0.051 


0.019 


0.17 
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0.542 


0.28 


0.41 


0.051 


0.018 


0.21 


2.2 


0.565 


0.39 


0.27 


0.048 


0.027 


0.22 


2.7 


0.589 


0.31 


0.38 


0.040 


0.012 


0.23 


3.05 


0.609 


0.50 


0.35 


0.003 


0.002 


0.10 


3.5 


0.664 


0.47 


0.33 


0.019 


0.013 


0.13 


3.75 


0.741 


0.48 


0.34 


0.0007 


0.0002 


0.14 



Specifically, the background of stellar models has been 
extracted from the evolutionary calculations recently pre- 
sented in Miller Bertolami & Althaus (2006) and Althaus 
et al. (2005), who computed the full evolution of initially 
1, 2.2, 2.7, 3.05, and 3.5 Mq. In addition, we include a new 
sequence of initially 3.75 M Q . All of the sequences were 
evolved from the ZAMS through the thermally pulsing 
and loss mass phases on the AGB. After experiencing sev- 
eral thermal pulses, the progenitors depart from the AGB 
and evolve towards high effective temperatures. Mass loss 
during the departure from the AGB has been arbitrar- 
ily fixed, as to obtain a final helium shell flash during 
the early white dwarf cooling phase. After the born-again 
episode, the hydrogen-deficient, quiescent helium-burning 
remnants evolve at constant luminosity to the domain 
of PG1159 stars with surface chemical composition rich 
in helium, carbon and oxygen. The masses of the rem- 
nants span the range 0.530 — 0.741M Q . For the sequence 
of IMq two different AGB evolution have been consid- 
ered, with different mass loss rates as to obtain different 
number of thermal pulses and, eventually, two different 
remnant masses of 0.530 and 0.542 M©. The main charac- 
teristics of the sequences considered in this work are given 
in Tabled We list the initial and final stellar mass (at the 
ZAMS and PG1159 stages, respectively), and the surface 
abundance of the main chemical constituents during the 
PG1159 stage. The sequence with initial mass of 2.7 M Q 
is the same presented by Althaus et al. (2005). 

Our PG1159 models are characterized by envelopes 
with uniform chemical compositions that extend from the 
surface downwards well below the driving region (i.e., the 
chemical composition at the driving region is the same 
than at the stellar surface). Thus, our models are not char- 
acterized by chemical gradients between the driving region 
and the stellar surface. Note that our post-VLTP models 
predict a range in the surface composition. In particular, 
the final surface abundance of helium spans the range 0.28- 
0.50 by mass 3 , which is in agreement with the range of 
observed helium abundance in most PG1159s (see Werner 

3 These abundances are not only determined by the stel- 
lar mass, but also by the deepness of the envelope convection 
and mass loss episodes after the VLTP (see Miller Bertolami 
& Althaus 2006 for comments). Also, the number of thermal 
pulses considered in the AGB stage and the efficiency of over- 



& Herwig 2006). Our sequences with helium abundances 
quite larger than the standard ones observed in PG1159 
stars will allow us to explore, at some extent, the role of 
helium in the instability properties of pulsating PG1159s. 
Mass loss episodes after the VLTP have not been consid- 
ered in the PG1159 evolutionary sequences we employed 
here. 

The pulsation stability analysis was performed with 
the help of a new finite-difference nonadiabatic pulsation 
code which is based on the adiabatic version described in 
Corsico & Althaus (2005, 2006). The nonadiabatic code 
solves the full sixth-order complex system of linearized 
equations and boundary conditions as given by Unno et 
al. (1989). Our code provides the dimensionless complex 
eigenvalue (lu) and cigcnfunctions (y±, • • • , ye) as given by 
Unno et al. (1989). Nonadiabatic pulsation periods and 
normalized growth rates are evaluated as II = 2n/^R.(a) 
and r] — — S(er)/5R(er), respectively. Here, 5ft(cr) and $s(cr) 
are the real and the imaginary part, respectively, of the 
complex eigenfrecuency a — (GM^/Rf) 1 / 2 w. Our code 
also computes the differential work function, dW(r)/dr, 
and the running work integral, W(r), as in Lee & Bradley 
(1993). In this work the "frozen-in convection" approx- 
imation was assumed because the flux carried by con- 
vection is usually negligible in PG1159 stars. Also, the 
e-mechanism for mode driving was neglected in the com- 
putations because nuclear-burning shells in PG1159 mod- 
els usually destabilize very short periods that are not ob- 
served in GW Vir stars (Kawaler et al. 1986; Kawaler 1988; 
Gautschy 1997). We employed about 3000 mesh-points to 
describe our background stellar models, most of them dis- 
tributed in the envelope region where all the pulsation 
driving and damping occur. We employed the "Ledoux 
modified" treatment to compute the Brunt- Vaisala fre- 
quency (N) (Tassoul et al. 1990). 

3. Stability calculations 

We analyze the stability properties of about 2400 stellar 
models covering a wide range of effective temperatures 
(5.4 ;> log(Toff) <; 4.8) and a range of stellar masses 
(0.530 < M*/M Q < 0.741). For each model we have re- 
stricted our study to unstable £ = 1,2 g-modes with pe- 
riods in the range 50 s <J n <J 7000 s, thus comfortably 
embracing the full period spectrum observed in GW Vir 
stars. In line with other stability studies of GW Vir stars, 
all unstable g-modes in our PG1159 models are driven by 
the K-mechanism associated with the opacity bump due 
to partial ionization of K-shell electrons of Cy-vi and 
Ovii-viii centered at logT w 6.2 (Quirion et al. 2004; 
GAS05). 

We start by discussing the stability properties of two 
template 0.589M Q PG1159 models. These properties are 
common to all PG1159 models of our complete set of evo- 
lutionary sequences. The normalized growth rate (77) in 

shooting play a role in the final PG1159 surface composition 
(Herwig 2000). 
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Fig. 1. The normalized growth rate in terms of period 
(in seconds) for overstable g-modes corresponding to two 
0.589M Q PG1159 models. The location of these models 
on the HR diagram is shown in Fig. [3] Circumscribed dot 
in panel (a) correspond to a mode with k = 70. 



terms of pulsation periods (II) for overstable I = 1 modes 
corresponding to the two selected models are shown in Fig. 

Model (a) is representative of the high-luminosity, low- 
gravity pre-white dwarf regime, and model (b) is typical 
of the low-luminosity, high-gravity phase, when the object 
has already entered their white dwarf cooling track (see 
Fig. O ■ Note that modes excited in model (a) have pulsa- 
tion periods in the range 1000 <J II <^ 2000 s, substantially 
longer than those excited in model (b) (100 < II <J 600 s). 
For each model, r\ reaches a maximum value in the vicinity 
of the long-period boundary of the instability domain. In 
other words, within a given band of unstable modes, the 
excitation is markedly stronger for modes characterized by 
long periods. This effect is particularly notorious in model 
(b), being the value of the growth rate for the shortest pe- 
riods more than seven order of magnitude smaller than for 
the modes with longer periods 4 . 

Fig. shows details of the driving/damping process 
in model (a) for two selected dipole modes. We restrict 
the figure to the envelope region of the model, where 
the main driving and damping occurs. Thick continuous 
curve corresponds to dW/dr for an unstable mode with 
k = 70, n = 1496 s and r) = 3.9 x 10~ 5 (marked as a 
circumscribed dot in panel (a) of Fig.^J, while the thick 
dashed one depicts the situation for a stable mode with 
k = 150, n = 3216 s, r\ = -1.23 x 10~ 3 . Also plotted is 
the run of the Rosseland opacity, k, and the run of its 
logarithmic derivatives, kt and k p /{T^ — 1). As can be 

4 We note, however, that for modes with the shortest periods 
the value of r] is so small (<C 10 -9 ) that they are only marginally 
unstable. 
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Fig. 2. The opacity and opacity derivatives (left scale), 
and the differential work function (right scale) for two 
selected £ = 1 modes in terms of the mass coordinate 
(lower scale) and the logarithm of the temperature (upper 
scale), corresponding to the PG1159 model (a) analyzed 
in Fig. □ dW/dr for the mode with k = 70 (k = 150) is 
normalized to an extremum value of +1 (—1). The vertical 
dot-dashed line shows the location of the maximum in the 
C/O opacity bump. 

seen, the region that destabilizes the k = 70 mode (where 
dW/dr > 0) is clearly associated with the bump in the 
opacity centered at — logg « 8.4 [q = (1 — M r /M*)], al- 
though the maximum driving for this mode comes from 
a slightly more internal region (—logg ks 8). Note also 
that in the driving region the quantity kt -I- n p / (T3 — 1) is 
increasing outward, in agreement with the well known nec- 
essary condition for mode excitation (Unno et al. 1989). 
Since the contributions to driving at — logg from 7.5 to 10 
largely overcome the damping effects at 6 < — logg <; 7.5, 
the mode with k = 70 is globally excited. At variance, the 
strong damping experienced by the mode with k = 150 
(denoted by negative values of dW/dr), makes this mode 
globally stable. The situation at the low-luminosity, high- 
gravity phase as in model (b) of Fig. ^ is very similar, the 
only important difference being that the driving/damping 
regions are located at considerably more external layers. 
This is due to an outward migration of the opacity profile, 
induced by the evolution of the star. 

3.1. The theoretical GW Vir instability strip 

Here, we examine the location of the unstable domains on 
the HR diagram. In Fig.|3|we show the evolutionary tracks 
for our complete set of PG1159 model sequences, where 
the thick portions of the curves correspond to models with 
dipole unstable modes. A well-definite instability domain, 
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Fig. 3. The HR diagram for our complete set of PG1159 
model sequences of 0.530, 0.542, 0.565, 0.589, 0.609, 0.664 
and 0.741M Q (from right to left). The portions of the 
tracks harbouring stellar models with overstable £ = 1 g- 
modes are emphasized with thick solid lines. Thin dashed 
(£ = 1) and dotted (£ = 2) lines are parameterizations 
of the blue and red edges of instability. The small white 
circles on the 0.589M© track show the loci of the two tem- 
plate models analyzed in Fig. ^ 



bounded by a red (cool) edge at high luminosities, and 
by a blue (hot) edge both at high and low luminosities, 
is apparent in the plot. The blue and red edges for dipole 
and quadrupole modes for each sequence are connected by 
thin curves as given by standard nonlinear least-squares 
algorithms. The instability domains for £ = 1 and £ = 2 
look very similar, although the edges for £ = 2 are slightly 
shifted to higher effective temperatures, and the region of 
instability is somewhat wider than for £ — 1. Fig.[3]should 
be compared with Fig. 5 of GAS05. The global agreement 
between our results and the predictions of GAS05 is ex- 
cellent, in particular for the sequence of O.589M — the 
only sequence in common between those authors and our 
work. 

At variance with GAS05, in this work we have em- 
ployed PG1159 models with different masses derived from 
the complete evolution of the progenitor stars. This has 
enabled us to extend the pulsational stability analysis to 
lower effective temperatures in the high-luminosity, low- 
gravity region. As a result, we have been able to found, 
for the first time, a reliable high-luminosity, low-gravity 
red edge of the GW Vir instability strip. Clearly, the red 
edge is markedly sensitive to the stellar mass, being more 
hotter for the more massive models. 

Our blue edge (both for dipole and quadrupole modes) 
cannot be exactly represented by a straight line. It is 




Fig. 4. Same as Fig.|21 but for the log(T e ff) — logg plane. 

clearly seen for dipole modes in Fig. |3J The departures 
from a simple linear relation have their origin in the differ- 
ent surface chemical compositions with which our models 
of different stellar masses reach the domain of the PG1159 
after emerging from the born-again episode. In fact, the 
unstable portions of the evolutionary tracks corresponding 
to models characterized by a surface helium abundance of 
« 0.28 — 0.39 by mass (that is, the sequences with masses 
of 0.542, 0.565 and 0.589M Q ) extend slightly beyond the 
linear parameterization of the blue edge, as compared with 
the case of the more massive models which have larger he- 
lium abundances in the envelope (0.50, 0.47 and 0.48 for 
stellar masses of 0.609, 0.664 and 0.741M Q , respectively) 5 . 
In fact, increasing the helium abundance at the driving 
region the efficiency of pulsational driving is reduced (see 
GAS05). This is in line also with the finding of Quirion et 
al. (2004) that decreasing the helium mass fraction at the 
driving regions, the blue edge of the instability domain 
shifts to higher effective temperatures. 

Since spectroscopic calibration of PGI159 stars gives 
effective temperatures and surface gravities, it is useful 
to see how the evolutionary tracks and the instability do- 
mains look on the log(T g) — \ogg plane. Fig. 01 plots the 
complete set of PG1159 sequences on that plane, in which 
we emphasize with thick solid lines the stages with over- 
stable £ = 1 g-modes. Again, we characterize the loci of 
the blue and red edges for dipole and quadrupole modes 
with thin curves as parameterized by standard nonlinear 
least-squares procedures. We postpone to 2] a complete 
discussion of the general agreement between our pulsation 
models and the observed GW Vir stars. 

5 Note that the helium abundance range predicted by our 
models (0.3 — 0.5) is coincident with the observed one for most 
of PG1159 stars, as reported by Werner & Herwig (2006). 
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Fig. 5. The instability domains on the log(T c ff ) — II plane 
for <?-modes with 1=1 corresponding to our complete set 
of PG1159 model sequences. Different line styles depict 
the situation for different stellar masses. 

In the current computations, the dipole unstable do- 
main for each evolutionary sequence is separated into two 
regions, one of them corresponding to the high-luminosity 
phase (low gravity) and the other corresponding to low 
luminosities (high gravity) (see Figs. |3| and |3J. The only 
exception is the sequence with 0.530M© which exhibits 
an unique unstable domain that extends uninterrupt- 
edly from the high-luminosity phase to the low-luminosity 
regime. This is at variance with the results reported by 
GAS05, who found that for all of their sequences (in- 
cluding that with O.53OM ) the tip of the evolutionary 
"knee" is pulsationally stable in the case of dipole modes. 
This discrepancy is connected with the different loci of 
the tracks on the HR diagram and on the log(T c ff) — logg 
plane. In fact, the evolutionary sequences considered in 
this work reach lower logg and luminosity values for a 
fixed log(T c ff) than in the ones of GAS05. This reflect the 
fact that in the present work we have considered PG1159 
evolutionary sequences derived from the full evolutionary 
computations of the progenitor stars. 

We explore now the ranges of periods of unstable 
modes. We begin by examining Fig. which shows the 
instability boundaries of dipole modes on the log(T e ff) —II 
plane for the complete set of evolutionary sequences. Note 
that generally the periods of unstable modes for each se- 
quence are clearly grouped into two separated regions, 
one of them characterized by long periods and high-radial 
overtones, corresponding to evolutionary stages before the 
knee at high luminosities, and the other one characterized 
by short periods (low k values) corresponding to the hot 
portion of the white dwarf cooling track (low luminosities). 
As mentioned before, the sequence with 0.530M Q exhibits 



an unique instability domain. Thus, this sequence shows 
instability even along the whole stages around the knee. 
The splitting of the instability domains into two separated 
regions can be understood in terms of the magnitude of 
the thermal timescale at the driving region. We refer the 
reader to the paper by GAS05 for an illuminating demon- 
stration of this. 

The high-luminosity domain of instability exhibits a 
strong mass dependence, being the longest unstable pe- 
riods shorter for the more massive models. In addition, 
for the less massive models the longest pulsation period 
is reached at lower effective temperatures as compared 
with the situation for the more massive models. In all the 
cases the long-period limit is attained shortly after the 
beginning of the instability domain. The shortest unsta- 
ble period for each sequence is also markedly sensitive to 
the stellar mass. Generally, the minimum period is smaller 
for the less massive models. Thus, the period-width of the 
instability domain is larger for the less massive models. 
As can be seen from Fig. the instability island on the 
log(T c ff) — II plane is almost vertical for the model with 
0.530M Q , but its slope gradually decreases as we go to 
sequences with higher masses. 

The low-luminosity domain, on the other hand, shows 
a moderate dependence on the stellar mass. The maxi- 
mum overstable period is always larger for the less massive 
models. The minimum overstable period, however, does 
not show a clear trend with the value of the stellar mass. 
We note that the shortest unstable periods are of sa 55 s 
and correspond to the sequence with 0.741M©. This short- 
period limit is substantially lower than that reported by 
GAS05, of about 190 s for their sequence of O.64M . It 
is important to note that in our computations the short- 
est overstable periods have very small growth rates, of 
the order of 10~ 9 — 1CP 13 , and thus we could consider 
that these modes are stable [see panel (b) of Fig. E]. If 
so, our short-period limit agrees with the predictions of 
GAS05. The morphology of the instability domains for 
quadrupole modes looks very similar to the case of dipole 
modes, being the novel feature a markedly shortening of 
the overstable periods, in agreement with GAS05. Indeed, 
the long-period limits of the high-luminosity instability 
domains for quadrupole modes are shortened in about 
3000 s as compared with the case of dipole ones. A less 
severe decreasing of overstable quadrupole periods is also 
present in the low- luminosity domain. As a result, both 
regions of instability are closer between them, to such a 
degree that the high- and low-luminosity instablility do- 
mains of the 0.530,0.542 and 0.565M Q sequences merge 
into an unique, uninterrupted region 6 . The lowest I = 2 
short-period limit for the low-luminosity domain, which 
corresponds to the sequence with M* = 0.741M Q , is of 
about 50 s. 



6 Note that in Figs. and |3| the blue edge for £ = 2 crosses 
the evolutionary track of O.565M0 ; we stress that it is simply a 
parameterization of the blue edge and not the exact theoretical 
blue edge emerging from our stability analysis. 
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In all of the computations presented in this paper we 
neglected chemical diffusion, and thus the stability calcu- 
lations were performed on PG1159 models with a constant 
chemical composition at the driving region. In particluar, 
no helium enrichment at the driving region was allowed, 
and consequently all of our sequences shown pulsational 
instability well beyond the empirical red edge of the GW 
Vir stars at low luminosities (T e g ~ 75 000 K; Dreizler & 
Heber 1998), even down to the domain of the variable DB 
white dwarfs. 

4. Theory vs observations 

In this section we compare our theoretical predictions with 
the observed properties of GW Vir stars. Currently, 11 
pulsating PG1159 stars are presently known. In Table|21we 
show the main spectroscopic and pulsation data available. 
Note that there are five GW Vir stars (termed PNNV) 
that are still embedded in a planetary nebula. The remain- 
der objects lack a surrounding nebula and are commonly 
called "naked" GW Vir stars. Finally, there are two ob- 
jects with measurable amounts of hydrogen in their spec- 
tra; they are termed pulsating "hybrid-PG1159" . Note 
that HE 1429-1209 is a naked GW Vir star but its ef- 
fective temperature and gravity place it at the region of 
the HR diagram usually populated by PNNVs. 

In Fig. we plot the location of pulsating and nonpul- 
sating PG1159 stars, as well as PG1159 stars that have 
not been observed for variability, on the log(T c ff) — \ogg 
plane (data taken from Werner & Herwig 2006). The plot 
also shows the evolutionary tracks for our complete set of 
sequences of PG1159 models. The blue and red edges for 
dipole and quadrupole modes are also shown. Regarding 
pulsating PG1159 stars, the agreement between observa- 
tions and model predictions is excellent. In fact, all of 
the GW Vir variables lay inside our predicted instabil- 
ity domains of dipole and quadrupole modes. We also see 
that, however, there are several non- variables occupying 
the unstable region. The existence of non-variable stars 
within the instability domain could in part be understood 
in terms of a variation in surface chemical composition 
(and thus in the driving region) from star to star. For 
instance, Quirion et al. (2004) found that the helium en- 
richment at the driving region is the cause for the exis- 
tence of the nonpulsator MCT 0130-1937 (with a helium 
abundance of about 75 %) within the instability strip. 
Note, however, that PG 1151-029, Longmore 3, Abell 21 
and VV47 (not included in the analysis of Quirion et al. 
2004) are found to have standard helium abundances (see 
Table 2 of Werner & Herwig 2006) and however are non- 
variables. On the other hand, it is remarkable the existence 
of the pulsating star NGC 246 with a helium abundance 
(Vhc ~ 0.63) unusually large among pulsators (see Table 
|2J). These controversial cases remain to be explained. 

Regarding PG 1159-035 (GW Vir) — the prototype 
of both the variable and the spectroscopic classes — our 
analysis naturally predict that this star does pulsate for a 
stellar mass of ~ 0.536M© (see Fig. EJ. It is a markedly 




Fig. 6. The distribution of the spectroscopically cali- 
brated variable and non-variable PG1159 stars in the 
log(T ff) — logg plane. PG1159 stars with no vari- 
ability data (N.V.D.) are depicted with hollow circles. 
Solid curves show the evolutionary tracks for differ- 
ent stellar masses: 0.530,0.542,0.565,0.589,0.609,0.664 
and 0.741M Q . Parameterizations of the theoretical dipole 
(dashed curves) and quadrupole (dotted curves) red and 
blue edges of the instability domain are also displayed. 



higher value than the mass required by GAS05 to excite 
pulsations in PG 1159-035 (M» w 0.50 - O.51Af ), and, 
at the same time, closer to the long-recognized "seismic" 
mass of w 0.59M Q of Kawaler & Bradley (1994). So, the 
employment of PG1159 models derived from the complete 
evolution of their progenitor stars appears to be a key 
factor to alleviate the discrepancy between the seismic and 
the spectroscopic mass of PG 1159-035. 

In what follows we shall focus on the observed pe- 
riod ranges for the known GW Vir stars. In Fig. we 
show the log(T c ff) — logn diagram, in which the effec- 
tive temperatures and the period ranges are taken from 
Table [5] For comparison, we have included the theoret- 
ical instability boundaries for 1=1. Also plotted are 
the ranges of unstable periods as predicted by Quirion 
et al. (2004). We stress that we are only interested in 
a qualitative comparison between the observed and the 
theoretical ranges of unstable periods. Thus, we do not 
attempt here a detailed asteroseismic period fitting for 
each individual GW Vir star. As can be seen from Fig. 
the agreement between theory and observations is rea- 
sonably good for all cases, and for several stars the agree- 
ment turns out to be excellent. Note, for instance, that for 
seven stars (Abell 43, HS 2324+3944, Kl-16, NGC 246, 
HE 1429-1209, PG 2131+066 and PG 0122+200) the ob- 
served ranges of periods are completely contained in the 
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Table 2. Stellar parameters and pulsation properties of all known pulsating PG1159 stars. The approximate surface 
abundances (in % by mass) have been derived from Table 2 of Werner & Herwig (2006) by assuming a composition 
made of 4 He, 12 C and 16 0, except for HS 2324+3944 (data taken from Table 1 of Werner & Herwig 2006) and for 
Abell 43 (data taken from Miksa et al. 2002). 
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theoretical instability domains corresponding to at least 
one stellar mass. For the case of PG 1159-035, the ob- 
served periods are well compatible with the instability 
domain of the 0.53Af Q models, although the shorter de- 
tected periods (<; 530 s) lay below the theoretical short- 
period limit. The opposite is true for PG 1707+427, whose 
longer observed periods (;> 760 s) lay above the highest 
long-period limit which corresponds, at the effective tem- 
perature of this star, to the 0.530M Q sequence. In the case 
of RXJ 2117+3412, the range of observed periods partially 
overlaps with four instability regions, corresponding to the 
sequences of 0.5890, 0.609, 0.664 and 0.741M Q , but no in- 
stability domain accounts for the observed periods shorter 
than w 800 s. Finally, in the case of Longmore 4 the band 
of observed periods is partially covered by the instability 
domains of the 0.530,0.542 and 0.741M Q sequences. Note 
that the agreement between the observed and predicted 
periods of unstable modes is comparable to that reported 
by Quirion et al. (2004). 

The predicted instability domains resulting from our 
analysis nearly reproduce the observational trend that 
the periods exhibited by pulsating PG1159 stars decrease 
with decreasing luminosity (increasing surface gravity) (see 
O'Brien 2000 for a discussion of such a trend). This is well 
documented by Figs. and Note that Longmore 4, be- 
ing a high-luminosity PNNV that shows short periods, is 
an exception to this trend. Longmore 4 is particularly in- 
teresting because it showed a surprising behaviour in its 
spectral type which suddenly changed from PG1159 to 
[WCE] and back again to PG1159. According to Werner 
et al. (1992), this could be a result of a transient but sig- 
nificant increase in the mass loss rate. Quite interestingly, 
according to our computations, Longmore 4 is located very 
close to the red edge of the instability strip (see Fig. EJ); 
so, this star could have just entered the instability phase. 
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Fig. 7. The period ranges of known variable PG1159 stars 
on the log(T c ff)— LT plane, depicted with thick vertical bars 
(the thickness of the bars is arbitrary). For comparison, we 
include the predicted ranges of unstable periods according 
to Quirion et al. (2004), shown with thinner vertical bars 
(arbitrarily shifted to right by 0.05 dex for clarity). The 
theoretical instability domains for (/-modes with I = 1 for 
our complete set of PG1159 sequences are depicted with 
different line styles. 

In addition to the issue of pulsation instability, valu- 
able information about pulsating PG1159 stars can be 
also extracted from the period spacing between consec- 
utive overtones of fixed I value. The period spacing de- 
pends primarily on the stellar mass and is only weakly 
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dependent on the luminosity and surface compositional 
stratification (Kawaler & Bradley f 994; Corsico & Althaus 
2006). Thus, this quantity allows a determination of M* to 
a very high accuracy. Here, we consider five GW Vir stars 
for which detailed asteroseismic studies have been car- 
ried out: PG1 159-035 (ATT = 21.5 s; Kawaler & Bradley 
1994), PG 2131+066 (ATT = 21.6 s; Reed et al. 2000), 
PG 1707+427 (An = 23.0 s; Kawaler et al. 2004), PG 
0122+200 (An = 21.1 s; O'Brien et al. 2000), and RX 
J2117+3412 (ATT = 21.5 s; Vauclair et al. 2002). We 
compare the observed mean period spacing of each star 
with the average of the computed period spacing, AIIj, 
(An^ = n^+i — life), corresponding to models with an 
effective temperature as close as possible to the value of 
T e ff of the star under consideration. 

We first consider the case of PG 1159-035. According 
to its location on the log(T e ff )— log g plane, this star should 
have a stellar mass of f» 0.536Af Q . We refer this value as 
the "spectroscopic mass" , M spoc (see Fig. |SJ . Our stability 
analysis predicts that a model with this mass at the effec- 
tive temperature of PG 1159-035 is pulsationally unstable. 
However, this model should have a value of An^ « 23 s, 
which is in conflict with the observed mean value of 21.5 
s. To have a An^ value compatible with the observed pe- 
riod spacing the stellar mass of PG 1159-035 should be 
of Msois ~ O.558M . Thus, we found a discrepance of 
AAI* = M sois — A/spec ~ 0.022M Q , somewhat lower than 
that reported in the literature (A/ spcc w O.54M , Dreizler 
& Heber 1998; M scis ~ 0.59M Q , Kawaler & Bradley 1994). 
We note that a model with M* « 0.558M Q at 140 000 K 
should be pulsationally stable in the frame of our stability 
analysis. 

In the case of RX J2117+3412, the evolutionary tracks 
of Fig.Elsuggest a stellar mass of about O.72M (see Miller 
Bcrtolami & Althaus 2006). A model with this mass and 
at the T e ff of this star should be pulsationally unstable. 
The corresponding value of An^ should be lower than 16 s, 
clearly at odds with the measured mean value of 21.5 s. So, 
in order to have An^ values comparable to the observed, 
we should be forced to consider models with masses lower 
than O.57M . Thus, for RX J2117+3412 we found a large 
disagreement between the spectroscopic and the seismic 
mass (AM* « — O.15Af ) and in the opposite direction 
than for PG 1159-035. Note that models with such low 
masses (at the effective temperature of RX J2117+3412) 
are outside of the instability domain (see Fig. EJ. 

Finally, we have the cases of PG 2131+066, PG 
1707+427 and PG 0122+200. According to our evolution- 
ary tracks, the mass of these stars should be of w O.55M 
for PG 2131+066 and of w O.53M for PG 1707+427 
and PG 0122+200. Note that, however, in order to have 
values of An^ compatible with the observed mean period 
spacings, the stellar masses should be substantially larger, 
of about O.58Af for PG 2131+066 (AM, w O.O3M ), 
0.57M Q for PG 1707+427 (AM* w 0.04M Q ), and of 
0.65M Q for PG 0122+200 (AM* w 0.12M©). Models with 
such high stellar masses (at the effective temperature of 
the stars under consideration) are pulsationally unstables 



(see Fig. |SJ). Given the spectroscopic uncertainties in the 
determination of log g, these solutions could be still com- 
patible with those derived from our evolutionary tracks 
and stability analysis. 

We conclude that the stellar masses of naked GW Vir 
stars as predicted by our evolutionary tracks are generally 
4 — 20% lower than those suggested by the period spac- 
ing data. By the contrary, for the PNNV RX J2117+3412 
the evolutionary tracks predict a spectroscopic mass about 
25% higher than the seismic derivation. Although our full 
evolutionary PG1159 models hint at generally agreement 
between the spectroscopic and seismic masses of pulsat- 
ing PG1159 stars, persisting discrepancies could still be 
reflecting a problem in the stellar modelling during the 
pulsing AGB phase of progenitor stars, as noted by Werner 
& Herwig (2006). 

5. Summary and conclusions 

In this paper we re-examined the pulsational stability 
properties of GW Vir stars. We performed extensive 
nonadiabatic computations on PG1159 evolutionary mod- 
els with stellar masses ranging from 0.530 to O.741Af . 
For each sequence of models, we computed the com- 
plete evolutionary stages of PG1159 progenitors starting 
from the Zero Age Main Sequence. Evolution was pur- 
sued through the thermally pulsing AGB and born-again 
(VLTP) phases to the domain of the PG1159 stars. The 
employment of such full evolutionary PG1159 models con- 
stitutes a substantial improvement over previous studies 
on GW Vir stars regarding the stellar modelling. 

Numerous detailed investigations about pulsating 
PG1159 stars have been performed on the basis of ar- 
tificial stellar models. In spite of the fact that the sig- 
nificant pulsation damping and driving occur in PG1159 
envelope stars, the employment of such simplified stellar 
configurations appear not well justified in the case of these 
stars. This is in contrast to the situation of their more 
evolved counterparts, the white dwarf stars, for which 
their thermo-mechanical structure has relaxed to the cor- 
rect one by the time the pulsational instability domains 
are reached. The main goal of the present work has been to 
assess to what degree the conclusions arrived at in previ- 
ous studies on PG1159 stars change when realistic stellar 
configurations are adopted. 

Our study confirms the following results, already 
known from previous studies: 

— g-modes in PG1159 models are excited by the n- 
mechanism due to partial ionization of carbon and 
oxygen. No abundance gradients between the driving 
region and the stellar surface are necessary to drive g- 
mode pulsations at the correct effective temperatures 
and period ranges. 

— There exists a well-defined instability domain with a 
blue edge which is strongly dependent on the stellar 
mass. 
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— Different surface helium abundances lead to sizeable 
differences in the precise location of the theoretical 
blue edge of the instability domain. 

— The instability domain splits into two separated re- 
gions, one of them at high luminosities characterized 
by long periods, and the other at low luminosities, 
corresponding to shorter periods, as recently demon- 
strated by GAS05. 

— All pulsating PG1159 stars lay into the predicted in- 
stability domain in the log(T e ff) — logg plane. 

— There is a very good agreement between the full period 
spectrum observed in GW Vir stars and the theoretical 
ranges of unstable periods. 

— The pulsation periods of excited modes decrease with 
decreasing luminosity (increasing surface gravity), in 
line with the observational trend. 

As for our new findings, we mention: 

— There exists a red edge of the instability domain at the 
high-luminosity (low-gravity) regime. This red edge is 
mass-dependent . 

— The border of the instability domains in the log T e g — 
log II plane at the high-luminosity, long-period regime 
is well delineated. 

— The pulsating PG1159 Longmore 4 is located at the 
very red edge of the instability strip at high lumi- 
nosities, a fact that could be reflecting the surprising 
behaviour observed in the spectral type of this star 
(Werner et al. 1992). 

— Some non-variables occupying the instability strip 
have standard helium abundances and the presence of 
them between pulsators can not be explained through 
the argument of Quirion et al. (2004). 

— The pulsating nature and also the range of observed 
periods of PG 1153-035 — the prototype of the GW 
Vir class — are naturally accounted for by pulsation- 
ally unstable PG1159 models with a stellar mass of 
- 0.53 - O.54M . 

Finally, we found that generally the seismic masses (as 
inferred from the period spacings) are somewhat different 
from the spectroscopic masses, although the disagreement 
for the PG 1159-035 star is somewhat alleviated accord- 
ing to our calculations. The persisting discrepancies could 
be attributed to a number of factors. On the observational 
side, possible systematics errors in the spectroscopic deter- 
mination of g and T e ff, and/or errors in the measurement 
of the period spacings of pulsating PG1159 stars. On the 
other hand, differences in the microphysics or the previous 
evolution may alter the location of the post-AGB tracks 
(Bloecker 1995). In fact, it has been argued by Werner 
k Herwig (2006) that the evolution during the TP-AGB 
(concerning third dredge up efficiency and TP-AGB life- 
times) may be key in determining the location of post- 
AGB tracks. However, in preliminary simulations we have 
found that neither third dredge up efficiency nor TP-AGB 
lifetimes play an important role in determining the loca- 
tion of post-AGB tracks. It remains to be seen if other 



physical assumptions like the overshooting efficiency dur- 
ing the core helium burning stage and early AGB (that 
also define the structure of the C-0 core and are com- 
pletely free parameters) may be playing a role in the lo- 
cation of post-AGB tracks. We are currently performing 
simulations of full stellar evolution sequences under differ- 
ent assumptions to clarify these issues. 
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